Relativistic Approach to Isoscalar Giant Resonances in 208 Pb 
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We calculate the longitudinal response of 208 Pb using a relativistic random-phase approximation 
to three different parameterizations of the Walecka model with scalar self-interactions. From a 
nonspectral calculation of the response — that automatically includes the mixing between positive- 
and negative-energy states — we extract the distribution of strength for the isoscalar monopole, 
dipole, and high-energy octupole resonances. We employ a consistent formalism that uses the same 
interaction in the calculation of the ground state as in the calculation of the response. As a result, 
the conservation of the vector current is strictly maintained throughout the calculation. Further, 
at small momentum transfers the spurious dipole strength — associated with the uniform translation 
of the center-of-mass — gets shifted to zero excitation energy and is cleanly separated from the 
sole remaining physical fragment located at an excitation energy of about 24 MeV; no additional 
dipole strength is observed. The best description of the collective modes is obtained using a "soft" 
parameterization having a compression modulus of if = 224 MeV. 



PACS numbers(s): 24.10. Jv, 21.10Re, 21.60.Jz 

Almost forty years ago Thouless wrote a seminal paper 
on vibrational states in nuclei in the random-phase ap- 
proximation jy . There he showed how spurious states — 
such as those associated with a uniform translation of the 
center-of-mass — separate out cleanly from the physical 
modes by having their strength shifted to zero excitation 
energy. Thirty years later Dawson and Furnstahl gen- 
eralized Thouless' result to the relativistic domain plac- 
ing particular emphasis on the role of consistency [Q. 
They showed how a fully self-consistent approach guar- 
antees the conservation of the vector current as well as 
the decoupling of the spurious component of the isoscalar 
dipole (J' r = 1~;T = 0) mode from the physical spec- 
trum. These results emerged after a careful treatment of 
the negative-energy states. Indeed, neglecting their con- 
tribution resulted in a violation of the vector current as 
well as the appearance of substantial spurious strength 
in the response. These fundamental results emphasize 
that the Dirac single-particle basis is complete only when 
positive- and negative-energy states are included. 

Relativistic models of nuclear structure have evolved 
considerably since they were first introduced by 
Walecka (U and later extended by Serot [Q. Although 
the qualitative success of these models relies almost ex- 
clusively on the dynamics generated by the the scalar 
(a) and vector (u>) mesons, several improvements have 
been introduced in order to enhance their quantitative 
standing |^,^| . Chief among them is the incorporation of 
scalar self-interactions which introduce important non- 
linearities into the equations of motion. Perhaps the 
greatest impact of these non-linear terms has been seen 
in the compression modulus of nuclear matter. In a lin- 
ear model the compression modulus is predicted to be 
unreasonably large at K — 547 MeV. Yet this value can 
be reduced to K = 224 MeV by the mere inclusion of 
non-linear terms. We will show here how this "soft" pa- 
rameterization yields excitation energies for various com- 



pressional modes in fair agreement with experiment. 

While scalar self-interactions are now incorporated 
routinely into most relativistic calculations of the nuclear 
ground state, their role on the dynamics of the excited 
states is just being unraveled; most calculations of the re- 
sponse of the mean-field ground state still use the linear 
model. Applying non-linear models becomes technically 
more difficult because the scalar-meson propagator no 
longer has a simple Yukawa form. At present the only 
calculations of the response that have incorporated non- 
linear terms are those by Ma and collaborators |?],|| . One 
of the main conclusions of their work is that "a large dis- 
crepancy remains between theory and experiment in the 
case of the dipole compression mode". We now show that 
if one includes the full momentum dependence of the lon- 
gitudinal response, a unique physical fragment emerges 
at low momentum transfer. This fragment — located at 
an excitation energy of E m 24 MeV — is identified as the 
isoscalar giant dipole resonance (ISGDR). 

We start from a Lagrangian having an isodoublet nu- 
cleon field (ip) interacting via the exchange of isoscalar 
sigma (4>) and omega (V 11 ) mesons, an isovector rho (b^) 
meson, and the photon (A^). That is, the interacting 
Lagrangian density becomes |||| 
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U(<j>) 



(1) 



In addition to meson-nucleon interactions the Lagrangian 
density includes scalar self-interactions of the form 



u(<f>) = + 1a0 4 



(2) 



Our theoretical program in the linear model has been 
described in great detail in several references j^Jl^] ■ Here 
we merely highlight the main features of the approach. 
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The longitudinal response of the mean-field ground state 
is defined by 
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where p(q) is the Fourier transform of the isoscalar vec- 
tor density, is the exact nuclear ground state, and 
is an excited state with excitation energy tv n . Note that 
the response is directly related to the timelike polariza- 
tion insertion IT 00 . To compute the linear response of 
the ground state of spherical nuclei — such as 208 Pb — one 
starts by calculating ground-state properties in a mean- 
field approximation. In this mean-field theory (MFT) 
nucleons interact with the self-consistent field generated 
by all positive-energy nucleons; vacuum loops are ne- 
glected in this approximation. Such a calculation yields 
single-particle energies and wave-functions for the occu- 
pied states as well as the mean- field potential Emf(i)- 
It is precisely this mean-field potential that one uses to 
compute the single- nucleon propagator nonspectrally: 

(w7°+i7 ■ V- M-Emf(x)) G f (x, y; u)=S(x - y). (4) 

There are several advantages in using a nonspectral rep- 
resentation for the nucleon propagator P, |ic||l2"f . First, 
one avoids the artificial cutoffs and truncations that 
plague the spectral approach ||. Second, both posi- 
tive and negative-energy continuua are treated exactly. 
As a result, the contributions from the negative-energy 
states to the response are included automatically. Fi- 
nally, a nonspectral evaluation of G-p poses no more 
challenges, nor requires much more computational ef- 
fort, than the corresponding calculation of an individ- 
ual single-particle state. Having determined the occu- 
pied bound-state orbitals and the nucleon propagator, 
the evaluation of the uncorrelated — or single-particle — 
polarization becomes relatively straightforward PJlCfl . 

To go beyond the simple single-particle response one 
must invoke the relativistic random-phase approxima- 
tion (RPA). The RPA builds long-range coherence among 
the many particle-hole excitations with the same quan- 
tum numbers by iterating the uncorrelated polarization 
to infinite order |pd| . Yet before going any further in 
the description of the RPA we must stress two issues of 
paramount importance. The first is consistency, which 
demands that the residual particle-hole interaction used 
in the RPA be identical to the interaction used to gener- 
ate the mean-field ground state. Second, the consistent 
relativistic response of the mean-field ground state in- 
volves, in addition to the familiar particle-hole excitation, 
the mixing of positive- and negative-energy states. These 
new configurations are essential for the conservation of 
the vector current and for the removal of spurious dipole 



strength from the physical spectrum. Although in the 
MFT it is consistent to neglect vacuum polarization [pj, 
the mixing between positive- and negative-energy states 
remains of utmost importance. 

The one new ingredient that we wish to add to our 
formalism is scalar self-interactions. The added compli- 
cation arises from the fact that the scalar-mediated inter- 
action no longer has a simple Yukawa form. Rather, the 
scalar propagator now satisfies a Klein-Gordon equation: 



(u 2 + V 2 -ml - U"(4>)) A(x, y; u) = c5(x - y) 



(5) 



In infinite nuclear matter the scalar self-interactions in- 
troduce a trivial modification: the scalar meson now 
propagates with an effective mass m* 2 = m 2 + U"(4>o), 
rather than with its free-space value. In the finite sys- 
tem solving for the scalar propagator becomes technically 
more difficult, but not more than solving for the nucleon 
propagator of Eq. (||). We have computed the scalar 
propagator in momentum space and have expanded it 
in terms of spherical harmonics so that the angular in- 
tegrals appearing in the RPA equations may be done 
analytically. A publication containing a more detailed 
description of our techniques will be forthcoming. 

The benchmark by which every theoretical calcula- 
tion of the nuclear response should be measured is the 
isoscalar giant dipole resonance. This is because the con- 
servation of the vector current and the shift of spurious 
strength to zero excitation energy can only happen in 
a consistent calculation of the response. In Fig. |l| we 
display the distribution of isoscalar dipole strength in 
208 Pb at the small momentum transfer of q — A6 MeV (or 
q = 0.23 fm _1 ) using parameter set NLC from Table 
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FIG. 1. Distribution of isoscalar dipole strength in 208 Pb at 
a momentum transfer of q = 46 MeV. Calculations were done 
using parameter set NLC |(| while the experimental value 
(filled circle on inset) is from Ref. [jwl. 

Note that the longitudinal response has been com- 
puted with an "artificial" width of 1 MeV. The un- 
correlated Hartree response displays a large amount of 
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dipole strength around 8 MeV of excitation energy. This 
strength is concentrated in the region where many 

particle-hole excitations can be made. Yet most of the 
strength is spurious, as revealed by the large amount be- 
ing shifted to zero excitation energy in the RPA response. 
What remains is an almost imperceptible fragment lo- 
cated at E = 24.4 MeV; and nothing else. The small 
fragment is displayed more clearly along with the experi- 
mental value (shown as a filled circle) on the inset of the 
figure. This result is a testimony to the power of con- 
sistency. By demanding that the residual particle-hole 
interaction be identical to the interaction in the ground 
state, and by properly including the mixing between 
positive- and negative-energy states, all spurious strength 
gets cleanly separated from the physical response. 

A comparison between three different relativistic 
models — all of them constrained to reproduce bulk prop- 
erties of nuclear matter at saturation as well as the 
root-mean-square charge radius of 40 Ca — is displayed 
in Fig. ||. Note that the three models employed here 
have been defined in Ref. j§ as L2 (K = 547 MeV), 
NLB (K = 421 MeV), and NLC {K = 224 MeV) [see 
also Table B. As expected, the energy of the dipole 
resonance scales with the compressibility of the model. 
Clearly, models with a large compression modulus — such 
as L2 and NLB — produce isoscalar dipole strength at 
values that are too large to be consistent with experi- 
ment fl3,[l4|]. These results have also been tabulated in 
Table H Because of the heroic efforts by experimen- 
talists in separating the isoscalar dipole mode from the 
high-energy octupole resonance (HEOR) , we also include 
a comparison between our results and their experimental 
findings in Table [n| Although not necessarily a com- 
pressional mode, our results for the HEOR follow similar 
trends to those observed for the giant dipole resonance. 
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FIG. 2. Distribution of isoscalar dipole strength in 208 Pb 
at a momentum transfer of q — 46 MeV for three different 
models. The experimental value is from Ref. pxj ]. 

We conclude the presentation of our results by dis- 



playing the distribution of strength for the quintessential 
compressional mode: the isoscalar giant monopole res- 
onance (GMR). First discovered in a-scattering experi- 
ments from 208 Pb p5| , and recently measured with higher 
accuracy at an excitation energy of E= 14.2±0.1 the 
GMR places important constraints on theoretical mod- 
els of nuclear matter. Indeed, the first measurement of 
the GMR — in conjunction with a a simple analysis based 
on the liquid-drop model — suggested a compression mod- 
ulus of about K = 200 MeV, a value considerably lower 
than the predictions of density-dependent Skyrme models 
at the time. Our calculations for the monopole strength 
in 208 Pb are displayed along with the experimental value 
in Fig. ||. We find good agreement with empirical for- 
mulas that suggest that the position of the GMR should 
scale as the square root of the compressibility. Indeed, we 
compute GMR energies in the ratio of 1:1.38:1.53, while 
the square root of the nuclear-matter compressibilities 
are in the ratio of 1:1.37:1.56. Moreover, these results 
help to reinforce our earlier claim that relativistic mod- 
els of nuclear structure having compression moduli well 
above K w 200 MeV will be in conflict with experiment. 
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FIG. 3. Distribution of isoscalar monopole strength in 
208 Pb at a momentum transfer of q — 129 MeV for three dif- 
ferent models. The experimental value is from Ref. |l6) . 

In summary, we have computed the distribution of 
strength for the isoscalar monopole, dipole, and high- 
energy octupole resonances in 208 Pb using a relativis- 
tic random-phase approximation to three different pa- 
rameterizations of the Walecka model with scalar self- 
interactions. We placed particular emphasis on the role 
of consistency. That is, we demanded that the residual 
particle-hole interaction used in the RPA be identical to 
the interaction used to generate the mean-field ground 
state. Moreover, we have used a nonspectral approach 
that automatically included the mixing between positive- 
and negative-energy states to compute the longitudinal 
response. Enforcing these constraints — and little else — 
was sufficient for separating the spurious J 11 = 1~ ; T = 
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state. In contrast to recent relativistic calculations ||— 
as well as nonrelativistic ones — we see no need 

for imposing additional constraints to "partially" remove 
the spurious contamination. These approaches attempt 
to remove all spurious strength by defining an effective 
dipole operator of the form Mio(r) = (r 3 — nr)Y w (r ); 
here 77 plays the role of a Lagrange multiplier and is 
determined to be 77 — 5(r 2 )/3 from translational invari- 
ance. More significantly, such a transition operator ne- 
glects the all important momentum dependence of the 
excitation. Indeed, it was only at small momentum 
transfers — just as in the experiment [|l4j — that we ob- 
served a single physical fragment concentrated around 
24 MeV of excitation. As the momentum transfer in- 
creased, we uncovered additional dipole strength around 
8 MeV. Yet this trend, namely, a sizable fraction of dipole 
strength at low energies and a giant resonance peak, is 
all that was reported in those recent publications flslJhSf. 
We do not regard this behavior as "a large discrepancy 
between theory and experiment". Rather, we attribute 
this trend to the simplified — momentum-independent — 
choice of dipole operator adopted in those calculations. 

We have also computed the distribution of strength for 
the giant monopole resonance. As in the case of ISGDR 
we have used the same exact operator — the isoscalar vec- 
tor density — to compute the monopole component of the 
longitudinal response. Indeed, monopole, dipole, and oc- 
tupole strength were all obtained from simply isolating 
the relevant J^-channel from the longitudinal response. 
For the GMR we have found good agreement with a re- 
cent relativistic calculation [Q] . Good agreement has also 
been obtained with semi-empirical formulas that suggest 
that the position of the GMR should scale as the square 
root of the compressibility. Depending on the relativistic 
parameterization adopted, monopole strength was found 
between 13 and 20 MeV of excitation. 

Lastly, we venture into the neutron-star domain. The 
non-linear parameterization NLC gives a rather satisfac- 
tory description of the various compressional modes. Al- 
though by no means perfect, this agreement suggests 
that the compression modulus of nuclear matter can 
not differ too much from the value predicted by this 
model (K = 224 MeV) . When this parameter set is used 
to compute the equation of state for neutron matter — 
and is then combined with the Tolman-Oppenheimer- 
Volkoff equation — one obtains an upper limit for the 
mass of a neutron star of M = 2.8M Q . Given the re- 
cent compilation of 21 neutron-star masses by Thorsett 
and Chakrabarty [^9) , where they show that the measure- 
ments are consistent with a remarkably narrow mass dis- 
tribution M = (1.35 ± 0.04) Af the fascinating pos- 
sibility that neutron stars harbor novel and exotic states 
of matter becomes almost a reality. 

This work was supported in part by the DOE under 
Contract No.DE-FG05-92ER40750 and by the Florida 
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TABLE I. Various relativistic parameter sets 
scalar mass and k are in MeV. 



The 



Set 


9l 


gl 


9l 


m s 




A 


L2 


109.63 


190.43 


65.23 


520 








NLB 


94.01 


158.48 


73.00 


510 


800 


10 


NLC 


95.11 


148.93 


74.99 


501 


5000 


-200 



TABLE II. Energies for various isoscalar resonances in 
three different relativistic models. All excitation energies are 
given in MeV. 



Model 



GMR 



ISGDR 



HEOR 



L2 
NLB 
NLC 



20.1 
18.1 
13.1 



31.3 
29.0 
24.4 



25.1 
23.4 
21.9 



Exp. 14.2 ±0.1 [16 22.4 ±0.5 [14 1 19.7 ± 0.5 [14] 
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